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PROBING SHOCK BREAKOUT AND PROGENITORS OF STRIPPED-ENVELOPE SUPERNOVAE 

THROUGH THEIR EARLY RADIO EMISSIONS 

Keiichi Maeda 1 

Draft 

ABSTRACT 

We study properties of early radio emission from stripped-envelope supernovae (those of type 
Ilb/Ib/Ic). We suggest there is a sub-class of stripped-envelope supernovae in their radio proper- 
ties, including optically well-studied type Ic supernovae (SNe Ic) 2002ap and 2007gr, showing a rapid 
rise to a radio peak within ~ 10 days reaching to a low luminosity (at least an order of magnitude 
fainter than a majority of SNe Ilb/Ib/Ic). They show a decline after the peak shallower than others 
while the spectral index is similar. We show that all these properties are naturally explained if the 
circumstellar material (CSM) density is low and therefore the forward shock is expanding into the 
CSM without deceleration. Since the forward shock velocity in this situation, as estimated from the 
radio properties, still records the maximum velocity of the SN ejecta following the shock breakout, 
observing these SNe in radio wavelengths provides new diagnostics on natures of the breakout and pro- 
genitor which otherwise requires a quite rapid follow-up in other wavelengths. The inferred post-shock 
breakout velocities of SNe Ic 2002ap and 2007gr are sub-relativistic, ~ 0.3c. These are higher than 
inferred for SN II 1987A, in line with suggested compact progenitors. However, these are lower than 
expected for a Wolf-Rayet (WR) progenitor. It may reflect a still-unresolved nature of the progenitors 
just before the explosion, and we suggest that the WR progenitor envelopes might have been inflated 
which could quickly reduce the maximum ejecta velocity from the initial shock breakout velocity. 

Subject headings: radiation mechanism: non-thermal - shock waves - circumstellar matter - super- 
novae: general - supernovae: individual (SNe 2002ap, 2007gr) 



1. INTRODUCTION 

Core-collapse supernovae (CC-SNe) are explosions of 
a massive star with the zero-age main-sequence mass 
(M ms ) exceeding ~ 8M Q . Their observational appear- 
ance is diverse, controlled by natures of the progeni- 
tor star, explosion, and the circumstellar environment. 
These are probably connected mutually through some 
intrinsic controlling factors, including at least the pro- 
genitor mass, metallicity, and the binarity. Clarifying 
these mutual relations is a key issue in the study of SNe. 
In this respect, stripped-envelope SNe (SNe Ilb/Ib/Ic) 
(Filippenko 1997) have been a target of various studies: 
Being an explosion of a star which has lost all or most 
of the hydrogen envelope before the explosion (Nomoto 
ct al. 1993), their origins are likely a mixture of differ- 
ent evolutionary paths (i.e., a single massive star with 
M ms > 25M Q or a binary evolution with M ms < 25M Q ). 
They show diverse properties in the explosion energies 
(Nomoto et al. 2010), with the most energetic ones linked 
to Gamma- Ray Bursts (GRBs) (e.g., Woosley & Bloom 
2006, for a review). In this paper, we address the follow- 
ing three issues in their natures. (1) The nature of shock 
breakout, (2) progenitor structure and its relation to the 
shock breakout signal, and (3) radio emissions and the 
CSM environment. 

The shock breakout emission is the first electromag- 
netic signal from SNe (Falk & Arnett 1977; Klein & 
Chevalier 1978). When the shock wave launched in the 
deepest part of the progenitor emerges from the stellar 
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surface, an intensive UV/X-ray flash comes out of the 
shock wave. The typical spectral energy of the emission 
is predicted to be sensitive to the progenitor radius, and 
for a Wolf-Rayet (WR) progenitor it is mostly in X-rays. 
The shock breakout phenomenon was predicted theoret- 
ically more than a few decades ago, but so far the detec- 
tion is rare (Soderberg et al. 2008; Schavinski et al. 2008; 
Gezari et al. 2008; Modjaz et al. 2009), first inferred from 
the ionization structure around SN 1987A (Lundqvist & 
Fransson 1996) - it is essentially a very brief transient 
event. This is now an important target of time-domain 
astronomy, including future optical survey proposals to 
catch the shock breakout signals from high-z SNe (e.g., 
Tominaga et al. 2011). However, such optical surveys 
are not optimized for SNe from a WR progenitor (i.e., a 
leading progenitor scenario for SNe Ib/c), since it is an 
X-ray event lasting only for ~ R*/c ~ 10 seconds (where 
i?„ is the progenitor radius, and c the speed of light). 

The final stage of the stellar evolution has not been 
clarified for stripped-envelope SNe. Estimating the ra- 
dius of a progenitor star provides important insight here, 
highlighted by the direct detection of an unprecedented 
blue supergiant (BSG) progenitor with i?„ ~ 50i? Q for 
SN II 1987A (see, Arnett et al. 1989, for a review). With- 
out the direct progenitor detection so far (see, Smartt 
2009, for a recent review), the shock breakout signal and 
the subsequent early optical emission have been used to 
estimate the size of the progenitors of SNe Ib/c (Soder- 
berg et al. 2008; Chevalier & Fransson 2008; Modjaz et 
al. 2009; Rabinak & Waxman 2011). SN lb 2008D is 
the only example where the shock breakout X-ray was 
convincingly detected (Soderberg et al. 2008; Modjaz et 
al. 2009), and the radius estimated here is still in de- 
bate (Chevalier & Fransson 2008; Rabinak & Waxman 



2 



2011). Any other independent information on the still- 
unresolved natures of the progenitors is highly valuable. 

Radio emissions have been detected from a number of 
nearby SNe Ilb/Ib/Ic (see, Chevalier & Fransson 2006b, 
for a review). The radio emission is created through the 
SN-CSM interaction, thus it is a strong tool to study the 
CSM environment. A relation between the time of the 
peak (t p ) and the luminosity there (L p ) can be used to 
estimate the size of the emitting region (thus the velocity 
of the shock wave at the SN-CSM interaction), and the 
CSM density (Chevalier 1998). Generally SNe of differ- 
ent spectral sub-classes occupy different regions in the 
tp — L p plot, and SNe Ib/c (plus a part of SNe lib, i.e., 
'compact' SNe lib: 'cllb' hereafter) are characterized by 
relatively low t p and large L p for given t p (Chevalier & 
Fransson 2006b; Chevalier & Soderberg 2010). However, 
radio emission from SNe Ilb/Ib/Ic is diverse in a sense 
that even excluding outliers they span more than an or- 
der of magnitude both in t p and L p (t p ~ 10 — 100 days, 
Lp ~ 10 26 - 10 28 erg s" 1 Hz" 1 ). Among SNe Ilb/Ib/Ic, 
SN Ic 2002ap is a suggested outlier characterized by small 
t p and L p . As we show in §2, SN Ic 2007gr shares similar 
properties with SN 2002ap, and these SNe show different 
features than others not only in t p and L p . 

In this paper, we propose a new idea to explain the pe- 
culiarities of radio emission from SNe Ib/c with low ra- 
dio luminosities and a short rising time (i.e., SN 2002ap, 
2007gr), and link our interpretation with the properties 
of shock breakout and progenitor. In §2, we summarize 
radio properties of SNe 2002ap and 2007gr, suggesting 
that these two belong to a sub class of stripped-envelope 
SNe (in terms of radio properties). In §3, we explore con- 
sequences of the hydrodynamic interaction between the 
SN ejecta and CSM in the early phase, and show that 
the evolution there can be different for SNe exploding in 
low density CSM environment (SNe 2002ap and 2007gr) 
and high density CSM (other SNe Ilb/Ib/c). In §4, we 
calculate multi-band radio light curves for SNe 2002ap 
and 2007gr, showing that their radio properties are nat- 
urally explained in terms of the hydrodynamic evolution 
discussed in §3. In §5, we discuss implications of our find- 
ings for natures of shock breakout and progenitor. The 
paper is closed in §6 with conclusions and discussion. 
Our scenario requires (relatively) inefficient acceleration 
of electrons at a shock, and discussion on this issue is 
given in Appendix. 

2. RADIO PROPERTIES OF SNE 2002AP AND 2007GR 

Properties of radio emissions from SNe are character- 
ized by a relation between the peak luminosity and the 
peak date (normalized at 5 GHz) (Chevalier 1998; Cheva- 
lier et al. 2006a). This is especially useful diagnostics 
for SNe Ib/c (Chevalier & Fransson' 2006b; Chevalier & 
Soderberg 2010): Since the synchrotron self- absorption 
(SSA) is a dominant absorption process for these SNe 
exploding within a relatively low density CSM environ- 
ment, a place in this t p — L p plot provides an estimate of 
the size of the emitting region (thus the velocity of the 
shock wave). The velocity of the shock wave estimated 
this way is typically ~ 0.1c for SNe cllb/Ib/Ic. At the 
same time, within the framework of the SSA, t p is smaller 
for less dense CSM. 

SN 2002ap is a classical example of broad-line SN Ic 
in optical wavelengths, suggested to be more energetic 



than canonical SNe (Mazzali et al. 2002). SN 2002ap 
is peculiar also in properties of radio emissions as com- 
pared to other SNe Ib/c (Berger et al. 2002; Bjornsson 
& Fransson 2004). While L p ~ 10 27 erg s^ 1 Hz^ 1 and 
£p(z/ p /5GHz) ~ 30 days for a majority of SNe Ib/c, in SN 
2002ap L p ~ 10 25 erg s' 1 Hz" 1 and i p (^ p /5GHz) - 3 
days characterized by an extremely low luminosity and 
a fast rise to the peak. According to the SSA scaling 
relation, the CSM density is estimated to be lower than 
other SNe Ib/c at least by an order of magnitude. In 
the optically thin phase, a majority of SNe Ib/c show 
the radio spectral slope a ~ — 1 and the temporal slope 
-!-5 < P < -1-3 (where L v oc v a t fi ) (Chevalier & Frans- 
son 2006b). However, SN 2002ap showed a shallow decay 
with f3 ~ —0.9 while the spectral index was similar to the 
others (a ~ —0.9) (Berger et al. 2002) (see also Bjornsson 
& Fransson 2004; Chevalier & Fransson 2006b). This pe- 
culiarity refuses a standard interpretation for SN 2002ap. 
There is only one theoretical interpretation suggested so 
far - a slope of the relativistic electrons' energy distribu- 
tion might have been different and flatter than in other 
SNe Ib/c (Bjornsson & Fransson 2004) (see also Cheva- 
lier & Fransson 2006b). 

We point out that SN Ic 2007gr shares the similar 
properties in radio with SN 2002ap, despite its optical 
properties belonging to a 'normal' (or non broad-line) 
class (Hunter et al. 2008; Valenti et al. 2008). In Figure 
1, we compare the multi-band radio light curves of SNe 
2002ap (Berger et al. 2002) and 2007gr (Soderberg et al. 
2010a). SN 2007gr showed L p ~ 10 26 erg s" 1 Hz" 1 and 
i p (f p /5GHz) ~ 5 days, both about an order of magni- 
tude smaller than typically found for SNe Ib/c, placing 
this SN close to SN 2002ap in these properties. Fur- 
thermore, the radio spectral index and decay slope in 
optically-thin phases are similar to those of SN 2002ap. 
Not only in t p and L p , but the temporal slopes arc differ- 
ent in SNe 2002ap and 2007gr (/3 > -1) from other SNe 
Ib/c (13 < —1.3), while the spectral indexes are similar 
(a~-l). 

The previous theoretical interpretation of the radio 
emission from SN 2002ap was based on the slow decay 
rate as mentioned above (Bjornsson & Fransson 2004), 
and this is our motivation to investigate an alternative 
explanation for this behavior. However, it should be 
noted that the quality of the radio data of SN 2002ap, 
as well as that of SN 2007gr, does not allow very accu- 
rate determination of the decay rate. Fitting the radio 
light curve of SN 2002ap (Berger et al. 2002) during 4 
- 50 days (8 points) by a function f v oc t@ , we obtain 
13 = -0.87 ± 0.17 at 8.46 GHz. 1 The error indicated 
here is lcr. For SN 2007gr (Soderberg et al. 2010a), we 
obtain /3 = -0.85±0.12 (4 - 18 days, 5 points). For com- 
parison, 'radio-normal' SN cllb 2011dh (Soderberg et al. 
2012; Krauss et al. 2012) shows j3 = -1.13 ±0.16 at 29.0 
GHz (30 - 100 days, 4 points). So, the preferred value 
of the decay rate (/3) is flatter by <~ 0.3 in SNe 2002ap 
and 2007gr than the canonical case. The deviation of 
SN 2002ap in the decay rate from the canonical value 
(J3 = —1.3 : Chevalier & Fransson 2006b) is therefore 
2.5a level, and for SN 2007gr it is above 3o\ Given the 

1 To ensure the optically thin nature of the emission, we compute 
the decay rate in the highest frequency band data available. 
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small number of the data points, however, this nominal 
significance should be regarded as merely indicative. In 
any case, in this paper we will investigate the implica- 
tions provided by this slow decay. 

As summarized in this section, SNe 2002ap and 2007gr 
(seem to) share the common properties, and we sug- 
gest they form a sub class of stripped-envelope SNe in 
the radio properties. Observationally this class of ob- 
jects are rare (about 10% of radio-detected stripped- 
envelope SNe), but they may well occupy intrinsically 
a large fraction of stripped-envelope SNe: SNe 2002ap 
and 2007gr are among the nearest examples of stripped- 
envelope SNe, and they would not have been detectable 
at ~ 30Mpc, i.e., the typical distance of most radio- 
detected SNe (Soderberg et al. 2010a). Future observa- 
tions of nearby SNe in radio wavelengths will be critical 
to establish if the shallow decay is a common property 
of these radio low-luminosity stripped-envelope SNe. 

3. HYDRODYNAMICS IN THE EARLY PHASE 

One of the striking features of radio emission from SNe 
2002ap and 2007gr is the short rise time with the time 
scale less than 10 days. In this section we discuss the 
hydrodynamic evolution of the shocked region in this 
early phase. In previous studies, the self-similar solu- 
tion for the interaction between the expanding SN ejecta 
and CSM has been assumed (Chevalier 1982). However, 
the basic assumption in the formalism would not apply in 
the very early phase of the expansion and/or low CSM 
density. The effect of the interaction on the hydrody- 
namic evolution becomes important only after a sufficient 
amount of CSM is swept up by the forward shock, and 
then the ejecta start being decelerated, following the self- 
similar solution which is eventually established. Before 
this phase, the ejecta feel the CSM almost as a vacuum, 
thus the ejecta are in the free expansion phase. 2 

The density structure at the outermost ejecta can be 
approximated by 



p S N~8.3 x 10 -18 
g cm" 3 



W Msn V 25 7 ^ V 101 Y 3 

(i) 



where psn is the density of the SN ejecta at the velocity 
V, E51 is the kinetic energy of the SN ejecta (in unit of 
10 51 erg), Msn is the ejecta mass, and td is the time since 
the explosion in day (Matzner & McKee 1999; Chevalier 
k Fransson 2008). ' 

The maximum velocity of the SN ejecta is determined 
by the shock breakout set by the radiative losses due to 
the shock breakout flash, unless there is any process that 
could alter the dynamics at the highest velocity ejecta 
after the breakout. Although detailed radiation hydro- 
dynamic modeling is required to obtain the exact value 
and the result depends on details in the treatment of the 
physics involved, the first-order estimate can be obtained 

2 The decelerated ejecta following the self-similar solution is fre- 
quently referred as the 'free expansion' in many literatures since 
the deceleration is not large in the self-similar solution considered 
here (i.e., R oc t 9 ) as compared to the Sedov phase (i.e., R oc t 0,4 
in CSM with constant density). In this paper, we refer the 'truly' 
free expansion phase (i.e., R oc t) before the establishment of the 
self-similar solution as the 'free expansion' phase. 



by analytic considerations (Matzner & McKee 1999): 
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where i?* is the radius of the progenitor (for that with a 
radiative envelope). 

Let us assume that the density distribution of the 
CSM is expressed by the steady wind solution, pcsm = 
M/4:Trr 2 v w (where M and v w are the mass loss rate and 
the wind velocity). Thus, 

Pcsm - 5 x 10 11 A*r~ 2 g cm" 3 , (3) 

where A* ~ 1 is a reference value corresponding to typi- 
cal Wolf-Rayet mass loss properties (i.e., M ~ 10~ 5 M Q 
yr -1 and v w <~ 1000 km s -1 ). Once the self-similar so- 
lution is established, the evolution of the velocity at the 
contact discontinuity follows the following form (Cheva- 
lier 1982; Chevalier & Fransson 2006b): 
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Since the self-similar solution describes the decelera- 
tion of the ejecta, the velocity here should be smaller 
than Vsn- Thus, the self-similar solution does not apply 
if 
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Assuming Vsn 3> 0.1c and the typical value of a few 
Mq for the ejecta mass and A* ~ 1, the strong interac- 
tion starts at latest a few days after the explosion, and 
thus the free expansion phase is negligible. This justifies 
the use of the self-similar solution for most of SNe. How- 
ever, if Vsn < 0.1c, then even a typical WR wind density 
is not enough to decelerate the ejecta to reach to the 
self-similar solution within the time scale of ~ 100 days. 
Furthermore, a low CSM density prevents the SN ejecta 
being decelerated, and with A* ~ 0.01 as inferred for 
sub class of SNe Ib/c (for example, the scaling relation 
of the SSA peaks implies that A, - 0.04 for SN 2002ap), 
the ejecta do not experience significant deceleration for 
~ 100 days even with Vsn ~ 0.2c. 

To confirm the above arguments, we have performed 
a series of hydrodynamic simulations for the early phase 
of the interaction. The ejecta structure is assumed to 
follow a power law up to the maximum ejecta velocity 
Vsn (equation 1), with E$i = 1 and M e j = 3M Q . The 
CSM density distribution is assumed to be a power law 
as described by equation 3. We have varied Vsn and 
^4* as parameters, choosing Vsn = 0.1c, 0.2c, 0.3c and 
^4* = 0.01,1,100 (thus 9 models were investigated in 
total). The interaction starts at 1 day since the explo- 
sion. Adiabatic Euler equations under spherical symme- 
try have been solved with the adiabatic index of 5/3, 
using the HLLC (Harten-Lax-van Leer Contact) solver 
(Toro 1999) to treat the discontinuity and the shock 
waves (Maeda et al. 2002). The number of meshes is 
8500 linearly covering up to 7.8 x 10 16 cm, so the spatial 
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resolution is about 10 13 cm which is sufficient to resolve 
detailed structure of the interaction region throughout 
the computation. 

Figures 2 and 3 show the evolution of the forward shock 
velocity and radius as a function of time. It is seen that 
the above analytical estimate roughly explains the be- 
havior of the shock wave velocity: Models with low Vsn 
and/or low A* do not follow the self-similar deceleration, 
but rather show an almost constant forward shock veloc- 
ity. According to equation 5, we expect idee as follows: 
For A» = 100, £dcc < 2 days for all values of Vsn, for 
-4* = 1, idoc < 1 day for V SN > 0.2c but ~ 200 days 
for V SN = 0.1c, and for A* = 0.01, t dcc ~ 2,60,20000 
days for Vsn = 0.3c, 0.2c, 0.1c, respectively. In the simu- 
lations, for A* = 100 the self-similar solution is approx- 
imately reached for all the cases. For A* = 1 the high 
velocity models with Vsn > 0.2c quickly follow the self- 
similar solution but the model with Vsn = 0.1c starts to 
be decelerated at t ~ 100 days. For A* = 0.01 the model 
with Vsn = 0.1c never experiences the deceleration dur- 
ing the simulated period of time, and the model with 
Vsn = 0.2c does not show a sign of deceleration before 
t ~ 40 days. 

Figure 4 shows the evolution of the velocity profiles 
for SN models expanding into the low density CSM with 
A* = 0.01, with the SN ejecta velocity V SN = 0.1c, 0.2c, 
and 0.3c. One immediately notices that all these three 
models do not follow the self-similar solution in the early 
phase, since in such a case there must be no difference 
according to different Vsn- In our simulations, the ve- 
locity of the leading edge of the ejecta is not identical to 
the initial value (i.e., Vsn)- At 6 days after the explosion 
(which may depend on the initial setup), the forward 
shock begins to develop, and the maximum velocity of 
the ejecta is - 75, 000 km s" 1 for V SN = 0.3c, 65, 000 km 
s" 1 for V SN = 0.2c, and 45,000 km s" 1 for V SN = 0.1c. 
Thus, the kinetic energy is redistributed in a way that it 
is decelerated for a large value of Vsn and accelerated for 
a small value of Vsn- This kinetic energy redistribution 
may partly stem from a specific simulation setup, and 
details of this process could be dependent on the initial 
conditions. Anyway, once the velocity profile is quickly 
adjusted, there is no further significant evolution here, 
and following arguments should not be sensitive to de- 
tails in the numerical setup. The forward shock is formed 
with the velocity roughly 80% of the (redistributed) max- 
imum velocity of the ejecta, and the forward shock ex- 
pands into the CSM at almost constant velocity without 
deceleration. 

It has been argued that the free-expansion phase is 
negligible for SNe Ib/c, with the shock breakout veloc- 
ity Vsn > 0.3c estimated with equation 2 or similar ex- 
pressions, under the assumption R* ~ Rq (Chevalier 
& Fransson 2006b). However, it has not been clarified 
yet to what extent the analytic expression on the shock 
breakout physics is accurate: the expression is an order- 
of-magnitude estimate with various assumptions. The 
progenitor structure just before the explosion is in de- 
bate, without direct detection of the SNe Ib/c progeni- 
tors so far. For example, a peculiar SN lb 2006jc showed 
a luminous blue variable (LBV)-like eruption about two 
years before the SN explosion (Pastorello et al. 2007). A 
mechanism of such an outburst of the proposed WCO 



Wolf-Rayet progenitor (Tominaga et al. 2008) has not 
been clarified yet, highlighting our ignorance of the final 
evolution of WR stars toward the SN explosion. 

Considering all these uncertainties in the shock break- 
out dynamics and progenitor scenarios, our motivation in 
this paper is to tackle to these issues from the opposite di- 
rection: We first provide possible, observationally-based 
information on these issues, then compare this with the 
theoretical expectations. Our strategy is the following. 
Rather than starting with the shock velocity based on 
equation 2 (or similar expressions), we first compare the 
velocity of the forward shock estimated from the radio 
data and that predicted by the self-similar decelerating 
solution. The two should agree if the dynamics is in the 
self-similar phase, while the radio-derived velocity will be 
lower than the self-similar expectation if the ejecta are 
in the free expansion phase. We then examine the radio 
spectrum and light curve properties of these SNe to see 
if the free-expansion solution provides a consistent view. 
Once we could obtain information on the free-expansion 
velocity through these analyses, then we could translate 
this information to the velocity of the shock breakout. 
This could then be used to calibrate the shock breakout 
physics (e.g., equation 2) and/or natures of progenitors. 

Useful diagnostics of the expansion of the forward 
shock, especially for stripped-envelope SNe, is provided 
by the observed peak date - peak radio luminosity re- 
lation. Figure 5 is a reproduction of the diagram from 
Sodcrberg et al. (2012) (see also Chevalier 1998; Cheva- 
lier et al. 2006a). Over plotted lines in the figure are the 
estimates which divide the (truly) free expansion phase 
and the self-similar deceleration phase using equation 5, 
for different choices of the ejecta mass. Below these lines, 
the observationally derived velocity (through the SSA 
model) is smaller than the self-similar prediction, thus 
SNe in this range are very likely in the free-expansion 
phase. 

Figure 5 shows that SNe Ib/c as well as cllb generally 
follow the self-similar solution with E$i ~ 1 and M e j ~ 
1 - 3M Q . An exception is SN 2002ap (as well as 1987A) 
which falls well below the self-similar prediction. Since 
different SNe have different properties in £51 and M e j, 
one has to compare the self-similar evolution and the 
radio properties on case by case basis. With E51 ~ 5 and 
M C j ~ 3M Q for SN 2002ap estimated through the optical 
modeling (Mazzali et al. 2002), the discrepancy between 
the self-similar solution and the observed radio properties 
is even larger. The same argument applies to SNe 2007gr 
(E 51 ~ 2, M cj ~ 2M Q ) (Hunter et al. 2008; Valenti et 
al. 2008) and 1987A (E 51 ~ 1.4, M ej - 14M ) (e.g., 
Blinnikov et al. 2000; Maeda et al. 2002), which showed 
that their radio properties are also below the self-similar 
case. Thus, we suggest that the shock wave evolution of 
SNe 2002ap and 2007gr (and 1987A) were in the free- 
expansion phase in these early epochs responsible for the 
early radio emission covering the SSA peak. 

4. RADIO EMISSION FROM SNE WITHIN LOW DENSITY 

CSM 

Results from previous studies on radio properties of 
SN 2002ap have been controversial. In the following, 
we show that a self-consistent picture can be obtained 
by considering the 'free expansion' which was not taken 
into account in the previous studies. We show further 
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that the same solution could apply to SN 2007gr as well. 
In our scenario, the peculiar features of SNe 2002ap and 
2007gr are natural consequences of SNe exploding within 
relatively low density CSM. 

We calculate radio emissions from SNe Ib/c as fol- 
lows. The basic formalisms have been developed by 
Fransson & Bjornsson (1998) and Bjornsson & Frans- 
son (2004), and specific prescriptions used here are given 
by Maeda (2012) (see also, Chevalier 1998; Soderberg et 
al. 2005; Chevalier et al. 2006a; Chevalier & Fransson 
2006b; Soderberg et al. 2010a). The synchrotron radio 
luminosity vL v in the optically thin phase is given as 



vL v ~ TTRl h V s hn rc ijl p m c c 



1 



synch(i^) 

t 



synch(i^) 
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(6) 

Here i? s h and I4h are the position and the velocity of 
the forward shock, n rc i is the number density of the rela- 
tivistic electrons. The relativistic electrons are assumed 
to follow a power law distribution with the index p as a 
function of the energy (note that throughout this paper 
we use p as the intrinsic power law index, before being 
altered by cooling effects). i sync h is the synchrotron cool- 
ing time scale, iothcr is the time scale for other energy 
loss processes which do not emit at the radio frequency 
(i.e., the IC scattering under the situation investigated 
here). The Lorentz factor of the electrons emitting at 
frequency v is 7„ ~ 80^°5 5 i?~ - 5 (here v w = v/10 10 Hz 
and B is in gauss). 

In SNe Ilb/Ib/Ic, the assumption of equipartition has 
been proven to work well (Fransson & Bjornsson 1998; 
Soderberg et al. 2005; Chevalier & Fransson 2006b). The 
energy densities of the relativistic electrons and the mag- 
netic field behind the shock wave are proportional to the 
thermal energy density created by the shock wave, while 
the proportional coefficients (e c and €b) are generally 
found to be lower than the full equipartition: Typically 
e c = €b = 0.1 is assumed for simplicity, while Fransson 
& Bjornsson (1998) and Maeda (2012) have suggested 
e c < 0.01 and es ~ 0.01 — 0.15. Our general arguments 
are independent from the values of e c and e^, for which 
further discussion is given in Appendix. The amplified 
magnetic field strength and the relativistic electron den- 
sity are expressed as follows: 



5-2.2 x lO^^A^-^gauss , 



, 0.5 _^sh_ 

Rsh 



n re i ~ 2.4 x 10 



17P-2 „ ( V sh 
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R 



cm 



(7) 
(8) 



Substituting these expressions into equation 6, the syn- 
chrotron emission is scaled as 



L v oc V sh v 



p/2 B (p-2)/2 



1 + 



^synch(z^) ^synch(i/) 



t 



t other {v) 



( 9 ) 

Note that so far no assumption has been made on the 
time-dependence of the expansion of the interaction re- 
gion. 

In these SNe, the main cooling agencies are the syn- 
chrotron cooling and the inverse Compton (IC) scatter- 
ing. These cooling time scales are estimated as 

W„chM ~ HOiv^B- 1 - 5 days , (10) 



ticH~1.7^- uo B 

Rsh 



0.5^0.5 ( L SN 



10 15 cm 



\ 10 42 ergs" 
2 

days . 



(11) 



Now, the spectral index (a) and the temporal slope (f3) 
of the synchrotron emission can be computed (here L„ oc 
v a t^). Describing the expansion dynamics as i? s h oc t m 
(therefore V s h oc i m_1 ), the result is summarized in Table 
1. 

A majority of SNe Ib/c show the spectral index a ~ — 1 
and the temporal slope —1.5 < (5 < —1.3 (Chevalier & 
Fransson 2006b) , being consistent with the adiabatic ex- 
pansion following the self-similar deceleration (m <~ 0.9) 
with p ~ 3. The 'low-density CSM' SNe Ic, 2002ap and 
2007gr, differ than the other cases in the temporal slope: 
they show a rather shallow decay, /? ~ — 1, while the 
spectral index is similar to the others (a ~ —1). One 
solution is to assume that the intrinsic slope of the rela- 
tivistic electron distribution is flatter (p ~ 2) than other 
SNe Ib/c (p ~ 3) and that the electrons emitting in the 
radio frequency are in the IC cooling regime (Bjornsson & 
Fransson 2004) , while assuming the self-similar dynamics 
(m ~ 0.9). In this case, the spectral slope is reproduced 
(see Tab. 1). The temporal evolution of the (optical) 
bolometric luminosity of SN 2002ap was roughly oc t° 
before 10 days since the explosion, and oc i~ ' 5 after that 
until 20 days (Yoshii et al. 2003; Tomita et al. 2006). 
Thus, if the IC cooling dominates, then the radio tempo- 
ral evolution will follow L u oc t 3 and oc t~ x before and 
after the optical peak (if the frequency of interest is op- 
tically thin) . Well after the optical peak, the cooling will 
become less important as time goes by, so the radio light 
curve will be flattened until it reaches to the temporal 
behavior of L v oc i~ - 8 (for p = 2). Thus, this scenario 
explains roughly the radio flux evolution, f3 ~ —0.9 as 
was observed. 

Although the scenario is consistent with the available 
observational data and also has a strength to explain 
the X-ray emission together with the radio (Bjornsson 

6 Fransson 2004), a drawback of the scenario is that 
it requires the electron distribution quite different from 
other SNe Ib/c. Note that for the magnetic field strength 
of B <~ 0.3 gauss (Bjornsson & Fransson 2004), the ra- 
dio emission is produced by electrons with the energy 

7 ~ 50 — 150, similar to those responsible for radio 
emission from other SNe cllb/Ib/Ic (e.g., Maeda 2012). 
Thus, the different slope cannot be attributed to differ- 
ent energy regimes probed for different objects. If this 
is true, it might mean that being a broad-line SN Ic, SN 
2002ap might have much more efficient electron acceler- 
ation (Chevalier & Fransson 2006b), but as we show in 
§2 the same argument should apply to SN 2007gr, which 
is a non-broad line canonical SN Ic, making this inter- 
pretation less appealing. 

In this paper, we suggest another solution for radio 
emission from SN Ic 2002ap, which also naturally ex- 
plains why 'canonical' SN Ic 2007gr shares the similar 
properties. The analysis in the previous section suggests 
that the ejecta expansion of SNe 2002ap and 2007gr is 
approximated by the free expansion without any decel- 
eration, rather than by the decelerating self-similar solu- 
tion. In this case, a ~ — 1 and (3 <~ — 1 are obtained with 
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p = 3, if the cooling is not efficient. This is the same 
situation as has been inferred for other SNe cllb/Ib/Ic 
(e.g., Chevalier & Fransson 2006b; Soderberg et al. 20f 2), 
with the only difference in the dynamic evolution: As- 
suming the general value of p = 3, then we predict that 
SNe in the low-density environment peak early in radio 
with the temporal index of B ~ — 1, while in the higher 
density SNe peak later and show the temporal evolution 
of B ~ -f.3. 

The radio light curves computed for the free expansion 
case as compared to those of SN 2002ap are shown in 
Figure 6. The microphysics parameters are set as eg = 
e e = 0.1 as typically assumed for radio SNe. As expected 
from Figure 5, we require relatively high forward shock 
velocity (V sh = 60,000 km s" 1 ) and low density CSM 
(A* = 0.007). The synchrotron cooling is included in the 
model, but its effect is negligible. Thus, the light curves 
show the temporal slope f3 ~ — 1 as observed (red-solid 
line in the figure), if the other cooling mechanism, i.e., 
the IC cooling, is ignored. 

Although it looks like a simple and straightforward in- 
terpretation, a story is complicated with the effect of the 
IC cooling. In Figure 6, we also show the light curves 
with the IC cooling included. With the set of parameters 
adopted, the IC cooling is indeed important, confirming 
the claim by the previous works (Bjdrnsson & Fransson 
2004; Chevalier & Fransson 2006b). With the standard 
p ~ 3, the model now fails to reproduce the spectral 
index. 

Since the match between the observed radio behaviors 
and the prediction without the IC cooling is striking, we 
investigate what conditions will be necessary to remedy 
the problem. The Compton cooling time scale is larger 
if e is smaller to reproduce a given luminosity. Figure 7 
shows an example where we adopt a low value for e e . The 
model parameters are V s h = 80,000 km s _1 , A» = 0.05, 
€b = 0.1, and e = 10~ 3 . The required mass loss param- 
eter is now increased roughly following the SSA scaling, 
e B A*(e c /e B ) 8 / 19 , but it is still very low as compared to 
other SNe Ib/c. The required velocity is also increased 
for the smaller value of e roughly following the SSA scal- 
ing 14h oc (e /es) -1 / 19 . Not only the decrease in the 
rclativistic electron density, but also the increase in the 
velocity, thus radius, has an effect to reduce the IC cool- 
ing effect. The light curve with low e is similar to the 
standard case. The difference is that in this model, the 
IC cooling is now negligible. With the parameters for 
this 'inefficient electron acceleration' model and E51 = 5 
and M C j = 3M Q , equation 5 predicts that tdec ~ 360 
days, justifying our assumption of the free expansion. 

For SN 2007gr, we have found the similar solution with 
SN 2002ap (Figures 8 & 9). The radio emission for the 
free expansion evolution explains the multi-band light 
curves fairly well without the IC cooling. With e c = 
€3 = 0.1 again the IC cooling alters the spectral index, 
and we require a low value of e e to fit the radio properties. 
For our fiducial model (V sh = 70, 000 km s" 1 , A* = 0.15, 
E 51 = 2, and M c j = 2M Q ), tdec ~ 40 days, which is long 
and roughly consistent with the observation. 

We note that the low value of e is not physically un- 
reasonable. Indeed, detailed modeling of radio proper- 
ties of SN ellb 1993J (Fransson & Bjornsson 1998) and 
SN cllb 2011dh (Maeda 2012) point to e c < 0.01, lower 



than generally assumed. No strong constraints have been 
placed for other SNe (e.g., Chevalier et al. 2006a). We 
discuss a few issues from the previous works on e e in Ap- 
pendix. The low efficiency of the electron acceleration 
has an important consequence in interpreting the X-ray 
emission, which is also discussed in Appendix. There, we 
conclude that so far there is no strong observational in- 
dication against the low value of e e in stripped-envelope 
SNe (and we suggest that the low value of e may be a 
generic feature). 

5. IMPLICATIONS FOR SHOCK BREAKOUT AND 
PROGENITORS 

5.1. SNe Ic 2002ap and 2007gr: Structure of WR Stars 
at The End of Their Lives 

Following our interpretation that the radio emission 
from SNe showing a short time scale to a low radio peak 
luminosity is described by the free-expansion of the SN 
ejecta without deceleration, we can potentially connect 
properties of the shock breakout and the progenitor to 
these early time radio properties. In this early epoch, if 
the forward shock is still in the free-expansion phase, the 
characteristic velocity obtained through the radio prop- 
erties is related to the maximum velocity of the SN ejecta 
following the shock breakout - thus, this important infor- 
mation can be obtained through the radio properties in 
the first month after the explosion, although in the other 
wavelengths it requires to catch the very moment of the 
shock breakout. We emphasize that this is not possi- 
ble for SNe already in the self-similar decelerating phase. 
Once the dynamics enters into the self-similar phase, the 
forward shock velocity is no more related to the initial 
maximum velocity of the ejecta (thus the shock breakout 
velocity) (see equation 4, also Fig. 2). 3 If there would be 
exactly the same two SNe except for the maximum veloc- 
ity (Vsn), discriminating these SNe requires information 

at t ^ tdec • 

According to §2 (Figures 2, 3, & 4), the freely expand- 
ing forward shock velocity is related to the maximum ve- 
locity of the SN ejecta after the shock breakout. Adopt- 
ing the parameters derived by the free-expansion model 
for the SN 2002ap radio emission, we have performed the 
same hydrodynamic simulation as in §2. We thereby find 
that V sh = 80, 000 km s" 1 is obtained if V SN ~ 90, 000 km 
s . In this simulation, the forward shock speed is con- 
stant until td ~ 40, thus being consistent with the free- 
expansion model for the radio emission. We have done 
the same experiment for SN 2007gr. Adopting A % = 0.15, 
indeed we find it is not possible to obtain the constant 
forward shock velocity of Vsh ~ 70,000 km s _1 for more 
than 10 days, thus this model (with e c = 10~ 3 ) for the 
radio emission is inconsistent with the dynamic evolu- 
tion. However, if we adopt A* = 0.05 (that corresponds 
to the required parameters of e ~ 0.01 and V s h ~ 60, 000 
km s _1 for ee = 0.1), then the constant forward shock 
velocity of V s h ~ 60,000 km s _1 is obtained for td ~ 50 
for the maximum ejecta velocity of Vsn ~ 75,000 km 
s _1 , being roughly consistent with the radio observation. 

Thus we infer the maximum ejecta velocity (or the 
'post-shock breakout velocity') of ~ 0.25 — 0.3c for SNe Ic 

3 It is much the same as the Sedov solution, where the informa- 
tion on the explosion is lost except for the explosion energy and 
the environment density. 
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2002ap and 2007gr through their early radio emissions. 
This is greater than inferred for SN 1987A, indicating 
that the progenitors of these SNe Ic are more compact 
than that of SN 1987A and/or the ejecta masses are 
smaller. However, the inferred velocity is smaller than 
the shock breakout velocity expected from the WR pro- 
genitors for SNe Ib/c (especially the WO or WC progeni- 
tor for SNe Ic). If ~ 5i? , one expects Vsn ~ 0.35c for 
typical SNe Ib/c, and ~ 0.6c for the ejecta properties of 
SN 2007gr, - 0.8 - 0.9c for SN 2002ap (adopting k ~ 0.2 
cm 2 g _1 for He or C+O composition). This could mean 
that there is still something missing in our understand- 
ing of the SN Ic progenitors and/or there is a missing 
part in our understanding of the shock breakout dynam- 
ics. In other words, the apparent discrepancy would pro- 
vide a hint on these still-unresolved issues, highlighting 
the importance of the independent information we could 
obtain through the early radio emission. We suggest the 
apparent discrepancy does provide an information on the 
progenitor structure just before the explosion. 

A possible solution we suggest is the following. While 
properties of typical Galactic WR stars are assumed 
in studying the SN explosion properties (including the 
shock breakout), we are not yet sure about what are 
properties of the WR star in a short period of time near 
the end of their lives. It has been suggested that the 
envelope would become expanded to a few to 10 times 
the original core radius, as a star approaches to the Ed- 
dington luminosity (Grafener et al. 2012). Studying SN 
properties so far has not provided any confirmation on 
this picture, and the new insight we obtain here could be 
the first indication that such an evolution could be the 
case. 

Grafener et al. (2012) showed that the typical density 
in such an envelope is ~ 10~ 10 g cm -3 , and ~ 1O _8 M 
is contained within ~ 10i? Q . The envelope can be more 
massive depending on the WR mass and other parame- 
ters, reaching to > 1O _6 -M in the model sequence stud- 
ied by Grafener et al. (2012). It is likely that this enve- 
lope would not affect the shock breakout itself - adopting 
£ 5 i = 1, M cj = 3Mq, and i?„ = 5 ft© , the stellar den- 
sity at the shock breakout is estimated to be ~ 10~ 8 
g cm~ 3 using the formalisms from Rabinak & Waxman 
(2011). 4 Thus, this would not dramatically affect basic 
predictions for the shock breakout high energy signal and 
subsequent early optical/UV emission, consistent with 
a few constraints favoring compact progenitors through 
the shock breakout X-ray (SN lb 2008D: Soderberg et 
al. 2008; Modjaz et al. 2009) and through the post-shock 
breakout optical emission (SN Ic PTFlOvgv: Corsi et 
al. 2012). 

However, the envelope could affect the dynamics just 
after the shock breakout. For typical parameters for SNe 
Ic, the shock breakout is estimated to take place when 
the mass above the shock (in the outer power law part, 
thus excluding the envelope here) is ~ 1O~ 8 M . This 
is indeed comparable to the envelope mass, or can be 
even smaller. Then, the highest-velocity ejecta will ex- 
perience deceleration during penetrating into the enve- 

4 However, the envelope could dilute and delay the shock break- 
out signal, as is similar to the situation expected for a shock break- 
out within a dense wind (Soderberg et al. 2008; Chevalier & Frans- 
son 2008). 



lope, before the velocity profile is frozen when the shock 
emerges out of the envelope. We suggest this is what 
we infer from the early radio emission. If the envelope 
mass is ~ 1O~ 6 M , it would decelerate the highest ve- 
locity to a half of the original shock breakout velocity, 
following Vsn/Vsn,o oc (IO^Mo/IO^M©)" - 15 (where 
^sn,o is the shock breakout velocity, while Vsn is the 
maximum velocity after being affected by the envelope). 
This relation was obtained by combining the equations 
presented by Rabinak & Waxman (2011). The envelope 
is probably required to be He-rich to make the envelope 
inflation, but this amount of He is much smaller than 
the upper limit obtained by spectral modeling of SNe Ic 
(Hachinger et al. 2011). 

5.2. Implications for Other classes of SNe 

The 'engine-driven' SNe linked to GRBs (Kulkarni et 
al. 1998; Soderberg et al. 2006a, 2010b) are distinguished 
from other stripped-envelope SNe in radio properties 
(Berger et al. 2003; Soderberg et al. 2006b). Accord- 
ing to Fig. 5, these radio-strong SNe are roughly consis- 
tent with arising from the forward shock wave following 
the self-similar expansion, with the SN highly energetic 
(E 5 i ~ 50). The observed luminosity is larger than the 
expectation nearly an order of magnitude, but we note 
that in this regime the relativistic treatment is neces- 
sary for quantitative comparison. The rough agreement 
suggests that the strong radio emission from these SNe 
may be understood in terms of the standard SN-CSM 
interaction scenario (with the relativistic ejecta which 
would require the 'central engine'), rather than invok- 
ing an additional relativistic jet component as long as 
the radio emission is concerned. The distinguished fea- 
ture in these engine-driven SNe is the large explosion en- 
ergy, which has been derived through the optical model- 
ing (e.g., Iwamoto et al. 1998; Woosley et al. 1999; Maeda 
et al. 2006). 

A majority of SNe Ib/c, and cllb are consistent with 
the self-similar expansion in the t p — L p plot. First, this 
indicates that their ejecta masses are mostly in the range 
of M e j = 1 — 3M Q . An important implication is that then 
a majority of them are likely an explosion of stars with 
M ms < 25M and require the binary companion to strip 
off their H-rich envelope: For example, if M ms ~ 25M Q 
(roughly a lower limit for a single massive progenitor 
for SNe Ib/c), then the ejecta mass is expected to be 
~ 4.5M Q if it explodes as an SN Ic and ~ 6.5M if an 
SN lb. Next, this sets a rough upper limit for the size of 
the progenitors for these SNe cllb/Ib/Ic. To enter into 
the self-similar phase, their shock breakout velocity must 
exceed at least 0.1c. For a reference value of E51 = 1 
and M e j = 3M , this means that the progenitor radius 
of most, if not all, of SNe cllb/Ib/Ic is < 25Oi? . 5 
Thus, we reject an RSG progenitor for SNe cllb/Ib/Ic - 
they must come from a WR star, or at most a yellow giant 
(the latter is a possible progenitor for SN cllb 2011dh: 
Maund et al. 2011; Van Dyk et al. 2011; Bersten et al. 
2012; Bcnvenuto et al. 2012). 

Another expectation is that some of SNe ellb ('ex- 
tended' SNe lib) may follow the free-expansion for a long 

5 The upper limit here would further go down if the structure of 
the progenitors of these SNe is also similar to what we suggest for 
SNe 2002ap and 2007gr. 
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time during the phases when the radio observations are 
performed, because of the expected low shock breakout 
velocity. With i?» ~ 500.R©, £51 ~ 1, and M cj ~ 3M , 
SN ellb 1993J is expected to have had the shock breakout 
velocity of 26,000 km s _1 . Taking this as a face value 
for SNe ellb, we expect that idee ~ 600,60,6 day for 
A* = 1, 10, 100. If there are SNe ellb within a relatively 
low CSM environment (e.g., A % = 10 for M = 10~ 6 M Q 
yr" 1 and v w = 10 km s _1 ), then such SNe ellb should 
show the free expansion phase in radio, and the high 
frequency observations will be especially useful to catch 
the feature. Indeed, SNe ellb tend to be below the 'self- 
similar expectation' in the t p — L p plot, inferring that 
some of them might be indeed in the free expansion 
phase. However, a complication is that in these cases the 
free- free absorption may become important, and either a 
more detailed model including the free-free absorption 
or the direct measurement of the shock velocity by the 
VLBI will be necessary. 

6. CONCLUSIONS AND DISCUSSION 

In this paper, we have investigated a consequence 
of SNe Ib/c exploding within low density CSM (i.e., 
A* < 0.1, or M < IQ- 6 M Q yr" 1 for a typical WR wind 
velocity). Such an SN should show a 'free expansion' 
phase before entering into the self-similar deceleration 
phase, if the post-shock breakout velocity is Vsn < 0.3c. 
The predicted radio properties are different from those 
SNe in the self-similar phase, characterized by a shal- 
low decline in the temporal evolution while the spectral 
index is the same, if all the other characters (i.e., the 
acceleration mechanism of electrons) are unchanged. 

SNe exploding in the low density CSM environment 
should be characterized in radio frequency by the fast 
rise to the peak (within 10 days after the explosion) 
and the low luminosity there. We have shown that all 
stripped-envelope SNe with these properties observed so 
far (SNe 2002ap and 2007gr, as well as SN 1987A from 
a blue giant progenitor) indeed have the expected prop- 
erties in the temporal and spectral indexes. The syn- 
thesized multi-band light curves show a good match to 
the observed ones for SNe 2002ap and 2007gr. We note 
that the other example, SN 1987A, was also well mod- 
eled by the free expansion dynamics (Chevalier 1998). 
In our interpretation, the efficiency of the acceleration 
of electrons must be low in order to avoid the IC cool- 
ing in radio frequencies. This is the same conclusion we 
obtained for SN 2011dh (Maeda 2012), and we suggest 
this may be a generic property shared by the SN-CSM 
interaction. 

Understanding the radio properties from SNe in the 
t p — L p plot in terms of the dynamics of the shock propa- 
gation, we propose new diagnostics on the shock breakout 
and progenitor properties through the early radio emis- 
sion: The forward shock velocity in the free expansion 
phase (V s h) is closely related to the maximum velocity 
obtained at/after the shock breakout (Vsn)- We sug- 
gest that the relatively low post-shock breakout velocity 
(~ 0.3c) we have derived for SNe Ic 2002ap and 2007gr 
could indicate the existence of an envelope driven by the 
near-Eddington luminosity near the end of the WR evo- 
lution (Graiener et al. 2012). This highlights the useful- 
ness/uniqueness of the proposed strategy as compared to 



other methods (e.g., breakout flash and/or subsequent 
optical emission) to probe the shock breakout and pro- 
genitor: the velocity information is not obtained by the 
other methods. Also, for a majority of SNe cllb/Ib/Ic, 
we reject a red-supergiant progenitor through the radio 
properties. 

The idea can in principal provide slightly different ap- 
proaches in estimating the post-shock breakout max- 
imum velocity and placing constraints on progenitor 
structures. Once one identifies the free expansion phase 
(i.e., the transition from the free-expansion to the self- 
similar dynamics in the decay slope) , then one could esti- 
mate the post-shock breakout velocity using equation (5) 
or similar expression, adopting ejecta properties and the 
CSM density estimated independently. The expression 
however requires further calibrations, and does not pro- 
vide a cross check of the decay slope. For these reasons, 
we adopted a more detailed approach, in which we fitted 
the multi-band light curves and checked the assumed dy- 
namic evolution with hydrodynamics simulations. Also, 
a stronger constraint on the progenitor radius when SNe 
do not show the free expansion phase than we did for 
most radio-detected SNe Ilb/Ib/Ic (using the informa- 
tion at the radio peak) could be placed, by using the 
earliest data points in which SNe are in the self-similar 
phase. A complication is that either the absorption (in 
the low frequency) or cooling (in the high frequency) can 
change the temporal slope in the early phase, thus distin- 
guishing the different dynamics is generally difficult well 
before the radio peak. For this reason, we have placed a 
conservative upper limit for the progenitor radius using 
the radio information around the peak. Further study 
and calibrations of these methods could be possible ap- 
plications of the idea presented in this paper. 

Future, large observational data set in radio frequency 
will be highly valuable. Such data are increasingly ac- 
cumulated recently thanks to great efforts of researchers 
working in the field (e.g., Soderberg et al. 2012). SNe in 
the low density CSM environment like SNe 2002ap and 
2007gr will provide a new possibility to tackle to proper- 
ties of the shock breakout and progenitor as mentioned 
above. A majority of radio-detected SNe Ilb/Ib/Ic follow 
the self-similar evolution, and thus their properties and 
distribution in the t p — L p plot will tell general distribu- 
tion of the progenitor mass and the energetics indepen- 
dently from optical wavelengths: So far, their distribu- 
tion suggests that a majority of them (i.e., radio-detected 
stripped-envelope SNe) are from relatively low mass pro- 
genitor (< 25Mq), indicating that a binary interaction 
path could be a dominant path, in line with other recent 
studies (e.g., Sana et al. 2012). 

The low density CSM around SNe Ic 2002ap and 
2007gr may be related to the WR progenitor structure 
just before the explosion. If the progenitor luminosity 
is larger, it would likely produce a higher velocity wind, 
resulting in a lower value of A* . This probably favors a 
massive WO star as a progenitor of these SNe Ic (e.g., 
Nugis & Lamers 2000). If this is true, we expect no 
SNe Ilb/Ib would belong to the 'rapid and faint' radio 
stripped-envelope SNe (i.e., SNe 2002ap and 2007gr are 
both of type Ic) , and future increasing samples should be 
able to address to this question. 6 Searching and observ- 

6 SN II 1987A may well be an exception, due to the low metal- 
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ing these radio faint stripped-envelope SNe will provide 
new clues on the progenitor scenario from this aspect as 
well. Especially, this will be best done by future observa- 
tories with the better sensitivity than currently. Because 
of the low radio luminosity of these objects, there may 
well be an observational bias in which we underestimate 
the frequency of these radio-weak SNe (Soderberg et al. 
2010a). Once a volume- limited sample is constructed, it 
will hopefully connect the radio properties with different 
progenitor evolutionally paths. 

licity environment there. This suggests that detailed analysis will 
require to take into account the metallicity of the local environment 
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Fig. 1 — Observed radio light curves of SNe 2002ap and 2007gr (Berger et al. 2002; Soderberg et al. 2010a). The behavior, L v oc u a t p as 
normalized at 1.43 GHz, is shown by lines (a = -0.85 and -0.9 for SNe 2002ap and 2007gr, respectively, and = -1.0 for both SNe). SNe 
2002ap and 2007gr show similar light curves, both in the spectral index and temporal slope, in the optically thin radio emission. Distances 
of 8 Mpc and 9.3 Mpc are adopted for SNe 2002ap and 2007gr, respectively. 




Fig. 2. — Evolution of the forward shock velocity computed by numerical hydrodynamic simulations, for VgN = 0.3c (green open squares), 
0.2c (red open circles), and 0.1c (blue open triangles). The CSM density parameter is A* = 0.01 (small symbols, top), A, = 1 (medium- 
sized symbols, middle), and A, = 100 (large symbols, bottom). The self-similar solution is shown by dotted lines for different values of A*. 
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Fig. 3. — Same as Figure 2, but for the evolution of the radius of the shock wave. See the caption of Figure 2 for meaning of different 
symbols. The self-similar solution is shown for different values of A* (dotted lines). The free-expansion evolution (i.e., R oc t) is shown by 
a solid line that fits well the evolution of the shock with Vsn = 0.2c and A* = 0.01 until ~ 70. 





Fig. 4. — Evolution of the velocity profiles at the beginning of the interaction, for the low density CSM with A* = 0.01. Three models are 
shown, with VgN = 0.3c (green), 0.2c (red), and 0.1c (blue). The position of the forward shock is indicated by the dashed lines, while the 
behavior expected from the decelerated self-similar expansion is shown by the dotted lines (with the vertical scale arbitrary). The model 
profiles arc shown starting at id = 1, with the interval of 2 days in different snap shots. 
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Fig. 5. — Radio properties of SNc Ilb/Ib/Ic and SN II 1987A in the t p — L p plot. The observed points are reproduction of Soderberg 
et al. (2012) (see also Chevalier et al. 2006a; Chevalier & Soderberg 2010), showing SNe Ic (cyan), lb (green), cllb (orange), ellb (red), 
'engine-driven' SNe Ic (including those associated with GRBs: blue), and 1987A (yellow). The velocity estimate with the SSA scaling is 
shown by dotted lines. The expectation from the self-similar hydrodynamic evolution is shown by the shaded area (blue, red, and green) 
for eg = 0.1 and a = e c /eg = 0.1 — 1, with E$i and M e j given in labels. The self-similar expectation for specific cases of SNe 2002ap 
(E S1 = 5, M oj = 3M ), 2007gr (E 5 i = 2, Af cj = 2M ), and 1987A (_E 5 i = 1.4, Af cj = 14) are shown by red-filled area. 
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Fig. 6. — Radio emission from SN 2002ap compared with the 'efficient electron acceleration' model. The model adopts V s h = 6.0 X fO 9 
cm s —1 , A* = 0.007, p = 2.7, e e = eg = O.f. The red-solid line is computed without the inverse Compton cooling, while the blue-dashed 
line includes the inverse Compton cooling with the bolometric light curve of SN 2002ap (Mazzali et al. 2002; Yoshii et al. 2003). The 
synchrotron cooling is included in both models, but its effect is negligible. 



(a) 2002ap, 1.43 GHz 






Fig. 7. — Radio emission from SN 2002ap compared with the 'inefficient electron acceleration' model. The model adopts V^h = 8-0 x fO 9 
cm s -1 , A* = 0.05, p = 2.7, t e = f0~ 3 , e B = O.f. See the caption of Figure 6. 




Fig. 9. — Radio emission from SN 2007gr compared with the 'inefficient electron acceleration' model. The model adopts V s h = 7.0 X 10 9 
cm s" 1 , A* = 0.15, p = 2.8, e e = 10~ 3 , eg = 0.1. See the caption of Figure 6. 
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TABLE 1 

Characteristics of the synchrotron emission 21 



Adiabatic Syn. IC 

a !-p _£ _P 

g (3m-3) + ^£ (3m -3)_+jg (5m - 5) + + <5 

"(P = 2) ^ =1 -1 

/3(p = 2) (3m - 3) - \ (3m - 3) (5m - 5) + 8 

<*(P = 3 ) ~ =1 =1 

/3(p = 3) (3m - 3) - 1 (3m - 3) - ± (5m - 5) - ± + 5 



a The spectral index (a) and temporal slope (/3) are shown for different 
cooling regimes (L u oc z; a f^). These are characterized by the electron 
distribution power law index (p) 1 the evolution of the forward shock (m, 
where i? s h oc t m ), and the evolution of the bolomctric luminosity (<5, where 
Lbol oc t 5 ). 

APPENDIX 

ELECTRON ACCELERATION EFFICIENCY (e E ) 

In our scenario, the acceleration of the relativistic electrons is required to be (relatively) inefficient (e < 0.01), as 
compared to frequently assumed (e ~ 0.1). As mentioned in the main text, so far there is no strong observational 
constraint against the low value of e . Indeed, Fransson & Bjornsson (1998) and Maeda (2012) suggested e < 0.01 
for SN ellb 1993 J and SN cllb 201 ldh, respectively, based on detailed modeling of their radio properties. We discuss 
this issue in this Appendix. 

Previous Works 

It has been sometime stated that e > 0.1 for SNe Ib/c based on the equation e e m p V^| ~ 7mra c c 2 where 7 m > 1 
is the characteristic Lorentz factor of relativistic electrons (Soderberg et al. 2005, 2010a). This is an analog of what 
is frequently adopted in the field of Gamma-Ray Burst non-thermal emission where the relativistic shock wave could 
accelerate all or most of the electrons into the relativistic speed. The assumption, that all the electrons are accelerated 
to the relativistic speed, however does not have to be the case in a non-relativistic shock wave in SNe. If only a 
fraction of Q (in number) are accelerated from the thermal population to relativistic energy, then one has to multiply 
the R.H.S of the above equation by Co, and therefore this limit on e e (by j m > 1) becomes lower by the same factor of 
( c . This is taken into account in our formalism (equation 8), and this is explicitly expressed as ( e ~ 4e ei _i(Vy0.1c) 2 
(from equations 3 and 8) for an H-rich CSM (of course it applies only for £ c < 1)- For reference values of e c = 10~ 3 
and V ~ 0.2c, we have Co ~ 0.1 — 0.2, i.e., about 10% of thermal electrons are accelerated to the relativistic speed. 

This value of ( c seems reasonable from circumstantial evidences: (1) £ e has been typically found to be less than 1% 
in SN remnants where V s h ~ 0.01c (e.g., Bamba et al. 2003). Thus we do not expect that the condition ( c ~ 1 must 
be met in the SN Ib/c shock wave which is also non-relativistic. (2) The X-ray emission from SN lib 1993J (similar to 
SNe Ib/c in the shock velocity) is believed to be emitted from the thermal electrons, and it requires that the thermal 
electrons are the bulk population at the shock front, while the relativistic electrons occupy only a small fraction in 
number (Fransson et al. 1996; Fransson & Bjornsson 1998). In any case, the value of ( e is still an open question (not 
only for SNe but also other astrophysical acceleration sites) . 

Furthermore, the existence of thermal electrons as a bulk population would not produce inconsistency in terms of 
observed radio properties. The synchrotron radio emission is produced by electrons with 7 > 50 for typical magnetic 
field strength seen in SNe Ib/c shock front, leaving no observational signature for electrons with 7 < 50 (including 
the thermal population) in the radio band (see Maeda 2012, for detailed discussion). Indeed, it is observationally 
forbidden to have j m » 1: In this case the radio emission should show a characteristic spectral break (Soderberg et 
al. 2005; Chevalier & Fransson 2006b; Soderberg et al. 2010a), although such a break is not seen in the observation. 
In this sense, the radio observation provides the upper limit for e e , not the lower limit. 

X-ray Production through the Inverse Compton Mechanism 

The inverse Compton (IC) up-scattering of SN optical photons has been proposed as one of the X-ray emission 
production mechanism, especially favored for stripped envelope SNe with A* ~ 1 (e.g., Bjornsson & Fransson 2004; 
Chevalier & Fransson 2006b; Soderberg et al. 2012). The scenario requires a large population of relativistic electrons 
with 7 < 50, and thus a large value of e (> 0.1) if these electrons' energy distribution follows a power law extrapolated 
from the radio-synchrotron emitting electrons (typically with 7 ~ 50 — 200) (see, Maeda 2012, for detailed discussion). 
In this section, we stress that having a large value of e (> 0.1) is not a single solution to produce the IC X-rays, but 
there is an alternative interpretation in which the energy distribution of the IC emitting electrons does not follow the 
extrapolation from the synchrotron emitting electrons as suggested by Maeda (2012). 

The low efficiency of the electron acceleration thus has an important consequence in interpreting the X-ray emission. 
Among several models suggested for the X-ray behavior of SN 2002ap, the IC scattering scenario seems the most 
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plausible (Bjornsson & Fransson 2004; Chevalier & Fransson 2006b). However, if we adopt p ~ 3, then we want to 
avoid the Compton cooling effect in the radio frequency (§4). We can estimate the IC X-ray luminosity in the form 
similar to the radio synchrotron emission: 



where t- lc is the IC cooling time scale. The typical electron energy (7; c ) where the Compton effect is significant is not 
sensitive to the value of e e since 7; c oc Lgj^-Rsh (not dependent on the microphysics parameters): At t<x = 5, 7; c ~ 140 
and 250 for our models with e c = 0.1 and 10 -3 , respectively. On the other hand, the electron energy responsible for 
the X-ray emission through the IC is 7 ~ 30. Thus it is in the adiabatic phase, and the predicted X-ray luminosity is 
v>L v (lkeV) ~ 6 x 10 36 erg s -1 (for e c = 0.1) and — 6 x 10 35 erg s _1 (for e c = 10~ 3 ) at day 5. The former is lower than 
the observed value only by a factor of two, but the latter is more than an order of magnitude lower than observed. 
The estimate is consistent with the previous work by Bjornsson & Fransson (2004). 

We note however that the similar situation is found in SN cllb 2011dh, for which the intensive radio and X-ray data 
allowed detailed modeling of its properties (Soderberg et al. 2012; Krauss et al. 2012; Bietenholz et al. 2012; Horesh 
et al. 2012). For SN 2011dh, there is almost no doubt that the power law index of the relativistic electrons emitting 
in the radio frequencies is p ~ 3, constrained by the late-time, optically thin light curves in multi bands (Soderberg 
et al. 2012). While Soderberg et al. (2012) suggested a high value of e e and Horesh et al. (2012) suggested even more 
efficient electron acceleration with e c /es ~ 1000, Maeda (2012) pointed out that little change in the spectral index was 
observed once the radio emission became optically thin, and suggested that the IC cooling effect must be negligible 
in radio and that the acceleration of the relativistic electrons must be inefficient (e c < 0.01). In this interpretation 
of Maeda (2012), the predicted IC emission in the X-ray was about one or two orders of magnitudes smaller than 
observed - this is the same situation we found for SNe 2002ap and 2007gr, but with the constraints on the spectral 
energy index of the electrons (p) and the effect of the IC cooling much stronger than for SNe 2002ap and 2007gr. To 
remedy this problem, Maeda (2012) suggested there is a distinct population of relativistic electrons below 7 ~ 50, 
with the number density (per 7) more than an order of magnitude larger than the extrapolation from the power law 
distribution for the radio-emitting, higher energy electrons, in order to explain the X-rays from SN 2011dh by the 
IC mechanism while still being consistent with the radio behaviors. We expect that the acceleration mechanisms is 
essentially the same in SN 2002ap, thus suggest that the same argument for SN 2011dh applies to SN 2002ap as well. 
Then the X-ray properties of SN 2002ap could be explained in the same manner. 
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